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_j Abstract 

— -^ In the article the problem of output setpoint tracking for affine non-linear sys- 

ps| tern is considered. Presented approach combines state feedback linearization and 

r^ homotopy numerical continuation in subspaces of phase space where feedback lin- 

^ earization fails. The method of numerical parameter continuation for solving systems 

of nonlinear equations is generalized to control affine non-linear dynamical systems. 

^^ Application of proposed method demonstrated on the speed and rotor magnetic flux 

control in the three-phase asynchronous motor. 



^ 1 Introduction 
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Let the affine nonlinear system with m inputs and m outputs in state space of dimen- 



cri sion n is given: 

> 

^ x = f{x) + ^gi{x)u,, y = h{x), (1) 

in *=i 

en where a; G X C M", y g F C M™, n G f/ C M™, maps / : M" ^ M", g, : M" -^ M", 

^ h : M" — i- M™ are smooth vector fields f,g,h G C°°. Functions f{x) and g(x) are 

considered as bounded on X. 

Systems of the form ([I]) are the most studied objects in the nonlinear control theory. 
k^ There are several most famous control methods for systems of type ([I| : feedback 

- linearization [I1I21E], application of differential smoothness [4j, Lyapunov functions and its 

generalizations [S] , including a backstepping [B] , also sliding control [7] and approximation 
of smooth dynamic systems by hybrid (switching) systems and hybrid control |HI. 

All of these control techniques have different strengths and weaknesses, their develop- 
ment is currently an active area of research, and the applicability and practical implemen- 
tation has been repeatedly con&med in laboratory tests and in commercial hardware. 

Approach described below is based on the method of numerical parameter continuation 
for solving systems of nonlinear equations [9], which deals with parametrized combination 
of the original problem, and some very simple one with a known solution. The immediate 
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motivation for the use of parameter continuation method in control problems is a series 
of papers [TOl [11], in which described the application of these methods directly in the 
process of physical experiments. 

In this paper we consider the solution of the output zeroing problem for the system 
([I| with relative degrees Vj > 1 that expands earlier obtained in [12] and [T3j results for a 
case Tj = 1. Further it is supposed that M it is free from zero-djTiamics, i.e. n = YlT=i ^j- 

The article consists of several parts. We briefly review the necessary facts about the 
method of parameter continuation and feedback linearization. Next, we represent the 
main result, an illustrative example of the method, as well as an example of controlling 
three-phase induction motor. 

2 Problem statement and motivation 

In this paper we consider the problem of nonlinear output regulation for affine non- 
linear system. In particular, we will solve the problem of output regulation to constant 
setpoint (without loss of generality, regulation to 0). 

Definition 1. Given the system of form ([I|. Problem, of output regulation to zero (aka 
output zeroing) is the design of such state-feedback control law u{t) = u{x) application of 
which asymptotically drives the system output to 0: limi_>oo?/(^) = 0. 

The output zeroing problem of affine nonlinear systems can be solved using mentioned 
above feedback linearization method. The main idea of the method consists in the trans- 
formation using a nonlinear feedback nonlinear system A^ : u{t) H- y{t) to the linear one 
L : v{t) I— )■ y{t) with the same outputs y, but new inputs v. After that, the resulting 
linear system L can be controlled by means of linear control theory. 

Suppose that a control problem of A^ can be in principle solved, i.e. there is exists 
a satisfying input signal u*{t), which gives the output response y*{t). The essence of 
problems in feedback linearization comes from that the response y*{t) may not be in any 
way reproduced by system L which is obtained after linearization. 

The simplest specific example is the system x = u, y = h{x) = x(a;^ — 1) + 1, a;(0) = 1 
for which the problem of output zeroing y -^ is needed to solve. 

If the system under consideration was a constant relative degree, the use of control 
V = —y after feedback linearization would give the output trajectory of y{t) = exp(— t), 
which is everywhere decreasing y{t) < 0. 

In this case, the nonlinearity y = h(x) has two limit points XI2 = il/v^, in which 
^a;(^i2) ~ 0- ^'^y trajectory y{t), that connects y{0) = 1 with y{T) = passes se- 
quentially through the points y^ = /i(3~^/^) and 1/2 = ^("3"^/^), and besides 7/2 > Vi- 
Hence, any trajectory y{t) on the interval (0,ti) should decrease with time (Figure 1), on 
the interval (^1,^2) increase, and in the interval again decrease. Such a trajectory is not 
reproducible using the feedback linearization. 
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Figure 1. Output trajectories of linearized system with constant relative degree and 

system with y = h{x) = x{x'^ — 1) + 1 

The behavior of the system in Figure 1 can be interpreted as follows: in the intervals 
(0,ti) and (t2,T) the system can be linearized in the usual manner and presented in the 
form y = V. On the interval (^1,^2) system behavior differs from the original, and the 
trajectory need to move in the opposite direction from the y = 0, which is the same as 
control of system y = —v. A similar situation arises in numerical methods for finding roots 
and optimization of functions with singularities, where in order to achieve optimum or find 
a root motion in the direction opposite to predicted by Newton's method is needed. We 
can use the parameter A G [0, 1] to indicate the motion direction. Increasing of parameter 
A > corresponds to the movement of y{t) in the direction to the desired setpoint y = 0, 
and parameter decreases A < in the opposite movement. The points of direction change 
A(t) = correspond to overcoming the singularities of h{x). In fact, this idea is the basis 
of the approach proposed below. 

3 Background 

In this section we present known facts needed to understand the main result. Finally, 
we come to the conclusion that the numerical homotopy methods can be used not only 
for solving nonlinear equations, but also for control of nonlinear afiine systems. Here and 
below we will always consider a setpoint tracking problem. 



3.1 Feedback linearization 

Definition 2. MIMO nonlinear system has relative degree Tj for output yj m © C M" if 
at least for one function gi is true 



c,^q-'h, ^ 



(2) 



where Cf\ = q^ fi^) ~ Yl'i=i dx- fii-'^) ^^ ^ -^^^ derivative of function A along a 
vector field f . 



It means that at least one input Uk influences to output yj after Tj integrations. 
Number r = YlT^i ^i i^ called as the total relative degree of system. If r = n and 
matrix 



/ Cg,CJ-^h,{x) ■ ■ ■ Cg^CJ-\i{x) \ 



A{x) 



(3) 



is full rank, then the original dynamical system ([l]) in © equivalent to system: 



2/f ^ = q h, + J2 l^sAl 'h, ■u, = B{x) + A{x) ■ u (4) 

The nonlinear feedback 



u = A{x)-^[v - B{x)\ (5) 

converts in subspace © original dynamical system ([T| to linear: 

y(^^) = v^ (6) 

Control of a nonlinear system ([T|^consists of two feedback loops, one of which imple- 
ments a linearizing transformation's]), second one controls the system ^ by any known 
method of linear control theory. 

A significant drawback, which limits the applicability of the feedback linearization in 
practice is requirement of relative degree r constancy and full-rank of matrix A[x) in the 
whole phase space &. 

3.2 Numerical continuation method 

Let it is necessary to solve system of the nonlinear equations 

0(0 = (7) 

where : W^ — )■ M™ is vector-valued smooth nonlinear function. 

Lets VL C W^ is open set and C{VL) is set of continuous maps from its closure Vt to W^. 
Functions Fq, Fi G (7(0) are homotopic (homotopy equivalent) if there exists a continuous 
mapping 

i7 : n X [0, 1] ^ M™ (8) 

that i/(^,0) = Fo(0 and H{^, 1) = Fi{Cl for all ^ G 0. It can be shown |9] that the 
equation H{^, A) = has solution (^, A) for all A G [0, 1]. The objective of all numerical 
continuation methods is tracing of implicitly defined function -ff(^. A) = for A G [0, 1]. 

Lets if : D — )■ W^ is C^-continuous function on an open set T> C ]R™+^, and the 
Jacobian matrix DH{^, A) is full-rank rank DH{^, X) = m for all (^, A) G P. Then, for all 
{^, A) G P exists a unique vector r G IR™^^ such as 

DH{^, A) ■ r = 0, ||r||2 = 1, det T^^JI' ^A > 0, (9) 

and mapping 

^ ■.V^R"'+\ ^ ■.{^,X)^T (10) 

is locally Lipschitz on V. 



Function (17) specifies the autonomous differential equation 



^ _M/(^,A), ^(0)=eo, A(0) = (li; 

which has a unique solution ^(t), A(t) according to a theorem of solution existence for 
the Cauchy problem. 

It can be shown jS] that integral curve 7(t) = (,^(t),A(t)) reaches the solution point 
(^*, 1) where 0(^*) = in finite time t < oo. 

In [T2| [T3] the method of numerical homotopy continuation for the time-dependent sys- 
tems of nonlinear equations 0(^, t) = is given. Another variant of numerial continuation 
for the system of nonstationary equations described in |15] . 

4 Main result 

In |12] pointed that the problem state control of the afiine system ([I| with relative 
degree of each state r^ = 1 associated with the solution ^*(t) : M+ — )■ M" for nonstationary 
nonlinear equation 0(^, t) = 0, in which g| = g{z), -^ = f{z). Namely, the control 

u{t) = J^ ' such that (f){^*{0),0) = a;(0), brings the system ([!]) from the state x(0) to 
x{T) = asymptotically T — )■ oo. Below we give another method than that described in 
that does not generate discontinuous control trajectories near limit points. 



4.1 Numerical continuation method for nonstationary system of 
nonlinear equations 

Consider the solution of the vector nonstationary equation (j){^, t) = 0, where : 
M" — 7- -R" is a smooth function. Compose parameterized simultaneously over time t and 
parameter A homotopy map: 

if (e, A, t) = (1 - A) ■ (^ - ^o) + A ■ 0(e, t) (12) 

where ^o is initial approximation to the solution. 

Lets formulate and give short proofs of some assertions for the background of parameter 
continuation method for nonstationary system of nonlinear equations. 



Assumption 1. If ^{t) and X{t) are solution of (12), then rank DH^^x{^{t), X{t),t) =m. 



From this assumption it follows that along the trajectory (^(t), A(t)) there is no bifur- 
cation points in which lank DH^x{^, X,t) < m and DxH G im D^H. In other words, the 
curve {^{t), A(t)) is free from branches and self-intersections. 

Equation H{^, A,t) = for A G [0, 1] defines implicitly defined function ^{t), parame- 



terized by A(t) and satisfies the equation obtained by differentiating (12) by time 



Denoting 



Bf-f rlT-f f)J-T 

Him, m,t) = -^m + -g^m + ^ = o (13) 



, 'dH dH\ ^ dH /f, 



(13) can be represented as a linear matrix equation with respect to ^ and A: 

At = B (15) 



Lemma 1. Equation (15) always has a solution. 



Proof Because oi A = DH^^x{^,^,t) the undetermined equation (15) has a solution if 



and only if rank DH^^x{^, A, t) = m, what is based on the assumption 1. D 



Lemma 2. All the solutions of (15) can be represented as t = a-r + r, where r = A'^B, 
T e kei A, a e R. 

Proof It is quite obvious, since dimkerA = 1 and all the null-space of A can be 
parametrized by scalar variable a G M, and the solution space W inhomogeneous equation 



of the form (15) is defined as H^ = {A+B} © kei A. D 



Let's prove assertion of Lipschitz maps \1' : I? — t- M^^^, which will need further. 

Theorem 1. Let if : "D — )■ M™ is C^ - continuous function on an open set T> C M™+^, 
and the Jacobian matrix A is full-rank rank A = m for all {C,, X,t) G V. Then for each 
(^, A,t) eV exists a unique vector f G ]R'^+^, such that 

f = a ■ T + f, f = A^B, 

f A\ (16) 

A-r = 0, ||r||2 = 1, det I ^ J > 0, a = const G M+, 

and map 

§:I?^M™+\ §: (e,A,t) ^f (17) 

is locally Lipschitz on T>. 

Proof Vector f = a ■ r + f defined as the sum of two components, one of which r is 
known (initial value problem 2.1.9 in j2]), that it is a Lipschitz function on T>. Hence, it 
is necessary to prove that f = A'^B = {DH^x)^DHtis Lipschitz. The uniqueness of the 
f follows from the uniqueness of the Moore-Penrose matrix pseudo-inversion. 

We can assume that DH is Lipschitz on V with a constant 7, which leads to existence 
and boundedness of the second derivatives H. Since rank A = m, the pseudo-inversion 
function A'^ continuous and differentiable [TT]. Hence the product A~^B is Lipschitz, 
because its components are Lipschitz. D 

Function \E' specifies the autonomous differential equation 



^ ,^, -^(e,A,t), e(0)=eo, A(0) = (18) 

which has a unique solution .^(t), A(t) according to a theorem of existence and unique- 
ness of solutions of the Cauchy problem. 

Theorem 2. The set if^^({0}) is simply connected. 



Proof By theorem 2.1 from ^0\, if for the map F : M" — > M'^, k < n true that 

sup{\\{DF{x)D'^F{x))-'^\\ I X G M"} < oo (19) 

then F~^{0) connected submanifold of dimension {n — k) b M". 

In accordance with the assumption [I| rank DH = m, therefore rank D^H = m, 
rank DH ■ D^H = m, inversion of the matrix {DH ■ D^H)^^ defined and its norm is 
bounded. D 

Theorem 3. A necessary condition for the existence of integral curve {C,{t), A(t)) for (18) 
that connects the points {^q, 0) and (^*, 1) is a > 0. 

Proof Based on lemma [2] and theorem [T] it is obvious that the Cauchy problem of the 
form (18) has a unique solution (^(t), A(t)), that satisfies the equation (12) for each fixed 
t = const: H{^{t), A(t), t) = 0. 

To determine the sign of a consider the behavior of the curve 7(t) = {^{t), A(t)) near 
t = 0. Since 

A=(A| + (l-A)£^-f).B=-A| (20) 



then near t = equation (12) behaves as a stationary with B = 0. Let's show that 
a ^ 0. If a = 0, then at time t = 0, from lemma ^^ and (20) it follows that ^(0) = 0, 



A(0) = 0, thereof A(t) = for each t and condition A = 1 never attainable. Since the 



near t = equation (12) stationary and the solution is determined by the conventional 



method of parameter continuation [9j, then orientation of vector r must be agreed with 
the condition I ^^ J > 0. Hence a > 0. D 

Now we state and give a short proof of the assumption about the behavior of the 



solution curve 7(t) = (^(t),A(t)) of (18). 



Theorem 4. There exists a number a^ G M., that the integral curve 7(t) = {C,{t), A(t)) for 



equation (18) with a > ao has finite length between the points (^OjO) and (^*, 1). 



Proof To prove this we consider the structure of the right side (18). Since the constant 



a can be selected arbitrarily large, then the term r can be neglected and may be written: 

f = a-T + f ¥a-T (21) 

By theorem [T] the function f is Lipschitz on V, then it follows automatically that f is 
bounded. Then one can always choose a finite a so that f = Q;-r + e, e^a-r. 



Since a is a finite number, then a ■ r as right side of ( 18 ) satisfies the well-known 
results on the finiteness of the solution trajectory (lemma 2.1.13 and theorem 2.1.14 in 
[6J). Hence, for an appropriate choice a > Uq curve (^(t), A(t)) has no limit points, and is 
diffeomorphic to the line, i.e. has finite length between the Aq = h Ai = 1. D 

The last theorem indicates that the parameter a is another one degree of freedom in 
designing the controller. The larger this constant, the faster the solution arrives to the 
A = 1, but numerical integration becomes more stiff. 



4.2 Homotopy continuation for nonlinear afRne systems 

Let's associate with the plant ([I]) hnear dynamics system with m inputs u, n states z, 
m outputs T] and with the same relative degries Tj for outputs such as in ([I]) 

z = Az + Bu, rj = Cz, -j^V = ^i- (22) 

Following equation is the homotopy mapping that links the outputs dynamics of the 
system (|l} and ( [22| ): 



i7= (l-A)-r/ + A-y = 0. (23) 

By definition of the relative degree of output, each component Hi should be differen- 
tiated Ti times with respect to t until it becomes an explicit function of any input u. We 
obtain after differentiation: 



r-i-l ri-1 



^"^ = - E ^'^^"'^AW + (1 - A)^, + (y. - r/,)A(-) + Y. C^^'-"^^^ 

"^ / . X '^' (24) 

+A f qhi + ^ Cg^Cf^K ■ ufc J = 0, 

that gives: 

i7(^') = A,i(a;,z,A) ■u + A,2(x,2,A) ■ A^^'^ + i3,(x,2, A), (25) 

where A = (A, A, A, ..., A*-^'"^^), C\ are binomial coefficients. 

Considering all of the components Hi after differentiation according to the relative de- 
grees of outputs Ti it is possible to write an algebraic condition that specifies a continuous 
deformation of system (22) to ([I|. 



H = Ai{x, z,A)-u + A2{x, z, A) ■ A^^'^^^ + I3{x, z, A) = 0, (26) 

where r^ax = max{ri}, A = (A, A, A, ..., A^'^™"-"^)), 

If i3 = 0, that is corresponds to the case r^ = 1, f{x) = 0, any known method of 
parameter continuation (like predictor-corrector method [9]) can determine the trajectory 

u{t), A(t) such as I • , , ) = r(t), where r(t) is tangent vector to the implicit curve H = 0, 

which is obtained from the linear matrix equation At = 0. This equation has infinitely 

many solutions as there are n conditions and n + 1 variables. In order to uniquely identify 

u{t) and A(t) an additional condition ||r|| = 1 for length normalization of the vector r 

needed. In addition, to select the correct direction of the r, imposed a condition of its 

fA\ 
positive orientation relative to the surface H, given in the form of inequality det I j^ J > 0. 

Considering the general case B ^ 0, tangent vector r(t) needs to be augmented by 
term for notstationarity compensating. Connected path {u{t), A(t)) for t G [0,T) starting 
at point {uq, 0) such as 

lim A(t) = 1, H{u{t), \{t),t) = 0, (27) 



can be generated from (26) as follows: 



u 



a ■ T + T, 



f = A^B, 



(28) 



A-T = 0, \\t\ 



1, det 



A 



>0, 



with the additional condition 



rank^(x, z,A) 



m. 



(29) 



where a > 0, a G M is a scalar constant, A'^ is Moore-Penrose inverse of matrix A. 



Condition (29) is a standard assumption when using parameter continuation method, 
which corresponds to the possible existence of limit points of trajectories ('u(t), A(t)) at 
which Ai ^ im^2; and the absence of bifurcation points. At the same time in some 
regions of phase space X x Z may be a situation where rank^i(x, z. A) < m, in that 
case, the system can not be linearized by the feedback, but the proposed method is 
applicable. Overcoming the bifurcation points, in which is observed Ai G im^2j also 
possible within the known approaches for the numerical parameter continuation (e.g., 
using the Lyapunov-Schmidt decomposition [9j). 



Equation (26) specifies the state feedback, and (28) specifies dynamics of the controller 



According to (24), equations for B and A2 depend explicitly on the y, and thus implements 



the output feedback. 



4.3 Switching strategies for nonregular feedback linearization 

All plants in practice are subject to a variation of parameters. The variant of control 
offered in the previous section is definitely sensitive to parametric uncertainties in a control 
object. On the other hand, in the feedback linearization control the parameter variation 
in the plant can be compensated by the controller for the linearized system |[16j. Let's 
consider a hybrid method that combines the possibility of applying an external control 
loop and resistant to change in the relative degree of the system. 



System in form (24) with output H can be linearized by feedback if consider evalution 
of parameter A as observable internal dynamics. In case if we fix A*^'''""'"-' = const, A*^'~'"°"='' G 



{— 1, 1}, then we obtain from (24) following affine nonlinear system: 



(30) 



H = F{x, z, A) + G{x, z,A) -u 

F{x, z, A) = A2ix, z, A) ■ A^*"™"^ + B{x, z, A) 

G{x,z,A) = Ai{x,z,A) 

It's possible to write a nonlinear coordinate transformation f i— )■ u transforming a 



nonlinear system (30) to linear one H^"^^ 



V 



u 



[Ai{x, z, A)]-i [v - A2{x, z, A) ■ A(^'— ) - B{x, z, A)) , (31) 

Switching strategy can be described as follows: we start whith conventional feedback 



linearization of (30), then in areas © where det (^(a;, 2;, A) f» and feedback linearization 



is not possible, it is necessary to switch control to the parameter continuation method 



implemented with (28) 



It should be noted that when v = control (31) is a special case of parameter con 



tinuation strategy that applied far from limit points. If A'-'"'"""''' 7^ 0, then control (31) is 



generated by equations (28) with following time-dependent scaling: 



sign(A('''""^)) — A^''™"^) f '^ \ f '"' \ - 

Since the degeneracy of matrix Ai in feedback linearization control always leads to 
increasing to infinity at least one input in u, practical way to perform switching between 
control strategies is consider event of inputs saturation when maxl-UjI > Umax- Input 
signals are always limited \u\ < Umax in the real world applications. This yields the 
following algorithm of the hybrid feedback linearization: we start with positive sign s = 1 
of A*^*'™"'"-* = s and conventional feedback linearization, then in case of saturation of inputs 
we switch to homotopy continuation procedure, and after return of all input signals to 
its limits we flipped sign s := —s and switch back to feedback linearization. Formal 
procedure for calculating the control actions can be represented as follows 

u := [Ai{x,z,A)]-^ (v - A2{x,z,A) ■ s- B{x,z,A)) 

if max \ui\ > Umnr then 



u \ (A 

),(._.)] =r, A-r = 0, ||r||2 = l, det(^^y 

u 



\(r ) I - A+B 

\ \i max } 

kx := (sign(A*^'''""^)) — A('"''^))/A*^'''"''^) 

ku := (sign(max \ui\)umax - max |Mj|)/max \ui\ 

if ku > kx then 

u := u ■ kx 

s := sign(A^''™''^)) 

elsefc„ > 1 

u := u ■ ku 

end if 

end if 

Since the set of regions & forms a compact space, using of this algorithm prodice 
almost linearized system. The situation where x E & can be considered as a perturbation 
action in the output H. 

5 Applications 

5.1 One illustrative example 

Consider following abstract example of MIMO system, that changes its relative degree 
in the state space 



10 



(33) 



Xi = Ui + X2 

yi = xl-xi + l 

1/2 = X2C0S(2X2) 

with initial conditions x(0) = (1, 1)^. We need to solve the problem of output zeroing 
y — > 0. Let's suppose that input signals constrained by inequality \ui\ < 20. 
Differentiating the outputs, we obtain 



yi 



{3x1 - 1) ■ (^^1 + ^2) = 9uUi + /i 



1/2 



(4^2 C0S(2X2) - 2^2 sin(2X2)) ■ {U2 + xf) = g22U2 + /2 



(34) 



Obviously, the system in interval x G [0, 1]^ can not be completely linearized by the 
feedback, because there are exists such x* that gn{x*) = or g22{x*) = 0. 



Let's associate with (33) linear system of a form 



Vi = ui, r]2 
with initial conditions 77(0) = (0, 0)"^ 



U2 



(35) 



According to the equation (|24|) we obtain for H = the following 

'gu 



Ai 



u + {l-X)E, A2 = y-v, ^ 



(36) 



7i 

The model in Simulink to control the system shown in figure 2. Modeling results are 
shown on figures 3-4. 
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Figure 2. Simulink model. 
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5.2 Three phase induction motor control 

Three-phase asynchronous motor is a famous example of a system that can not be 
linearized by state feedback |17| . Consider the application of the proposed method to 
control the speed and flux linkage of the motor. For modeling of the electric motor in the 
state space, we strictly follow the material of the paper [18j. 

Let's consider the reduced fourth-order state-space model of induction modor: 



12 



U = p 






-r.„ ^(hrd + r.„ ^M^ri 






«,s 



-/3a;0r - n ^isq - ^s^sd + 



"sq = -PUJ(pr - T^ ^Isq " tOs^sd + y^ 

J-il 

with this kind of parainetrization: 



Tr= -f^, li = P -f—, P 






Li — Lg -^, Ri — Rg + Rr [ — — > Ti „ 



(38) 



where isd and i^g are respectively the stator currents projections on the (d, q) axis 
reference frame, 0,. is a rotor fluxe, L^ and Lr are the stator and rotor self-inductances 
and Mgr is the mutual inductance. 

The electromagnetic torque developed by the motor is expressed in terms of rotor 
fluxes and stator currents as: 

Tern = P^^i'isq4>r) (39) 

where p is a number of pole pairs. 

Synchronous rotor angular speed Ug can be expressed as 

uJs = UJ^ -— (40) 

TrCPr 

The outputs to be controlled are the mechanical speed yi = u/p and the square of 
the rotor flux magnitude 1/2 = (pl- State variables are stator currents {isd,isq), the rotor 
fluxes {(prd, 4'rq) cLud the rotor angular speed u. Control variables are stator voltages Vgd 
and Vsq. 



First differentiation of outputs yields 

ijl = u/p = fKprisq/P - Tm/J 

After second differentiation of outputs finally inputs appeared: 



yi\ ^ (bi\ ^f ai2 



(41) 



61 = -fl(j)r{isq'Ti ^ + isqTr ^ + ^s^sd) + /^^r ^Mgrisdisq " /i/3w0r 

^(2 + /3M,,)0^ - (3^ + ^) ^^^ + 2MgrCOs^^ + 2Ml^ (42) 

T^ \ Tr Ti J Tr T^ T^ 



Q'12 — l^i^l 
^21 



2Msr<Pr 



TrLi 



Plant model in the form of (42) can be linearized by feedback when 0^ 7^ 0, after 



nonlinear transformation we will have 
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y = v (43) 

In [18] to control the (43) proportional-differential (PD) controller used, which in 
practice has a number of fundamental problems of reducing the stability to noise in the 
feedback. 



In this paper we propose a different approach to the control of (42), based on two feed- 



back loops: the internal to stabilize the current isd, isq and the external to control outputs 
yi,y2- As a result, only proportional-integral (PI) controllers are used and structure of 
the system resembles a classical FOC-control with the only difference being that the out- 
put of each PI controller is passed through an appropriate nonlinear transformation of 
coordinates. 

Parameter continuation is used only in the outer control loop, the inner loop is imple- 
mented with a current isd, isq decoupling by coordinate transformation |21) . 

Inner loop for current stabilization is implemented using a nonlinear feedback through 
which control signals z/^^, Ugq are passed 



Usd = LAusd -(^si 



sqj 



U 



sq 



Lliiysq + P(^<Pr +(^sisd) 



(44) 



which gives the decoupled linear dynamics of the currents 



^sd + /3l"r Vr - T-i ^i 



sd 



(45) 



"sq 



'sq 



' 1 ''sq 



Control of (45) can be achieved with a simple PI controller 



V = Kpif^f -i) + K, / ii'^f^r) - i{T))dT 



(46) 



where i^'^^ is a current setpoint for the corresponding axis. 

With corresponding adjustment of coefficients Kp and Ki can be achieved fast reg- 
ulation and exact match 
current regulation |21]. 

Let us turn to the outer loop to control the outputs of yi = (pl and the mechanical 



■ref 



~ ''sd' , 'isq ~ ilq'^ , which allows ignore the dynamics of the 



speed ^2 = ^/p, whose dynamics is given by the equation (41). Let's associate with (41) 
linear system of a form 



r]i = ui = isd, r]2 = U2 



''sq 



(47) 



with initial conditions 77(0) = (0, 0)"^. 



and (47) 



Using (24) and (25) we can write the following equation for mixed dynamics of (41) 



H 



Hi 



A(-2r-V' + 2M,,r-i( 







(48) 



Xij + n(l -\) + \{y-vi) = Q 

r-isd) + isd{l - A) + \{yi - r]i] 

H2 = Xfi4>risq/P + isq{l - A) + A(y2 - ^2) = 

The purpose control is the asymptotical output zeroing Hi{t) — )■ 0, H2{t) — )■ 0. Should 



be noted that we droped therm T^/ J in equation (41) for ^2- J\^^ SToro ocymecTBHM 



jiHHeapH3an,Hio chctcml.! (48) 
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Equation (48) establishes algebraic condition for the continuous deformation of the 



system (47) to (41 ). It is used to control the system in regions where feedback linearization 



is not possible, i.e. if 
form of 



7^ 0. The control inputs calculated by (28), which has particular 



I'sd 



2\(f)f. + ViTr + r]\Tr — XTrlJ 
Tr — XTr + 2MsrX(f)r 

V2 + r]X- Xy 



''sq 



flX(f)r/p + 1 — A 

with initial dynamic of paremeter A: A(0) = 0, A(0) 
controlled with convential PI regulators 



(49) 



1. Inputs vi and V2 are 



ref 



(50) 



fi = i^pid/^'^ - yi) + Kii ■ d, ei = yf - yi 

V2 = Kp2{yl''^ - y2) + Ki2 ■ Ca, 62 = 2/2"^ - y2 
When the system is far from area of singularity of linearizing transformation (i.e. when 



> 0) control (49) reducec into (51) 



Tr 



2M 
i 



ivi + 2r, 



-1 J,2\ 



sq 



P 



(51) 



-V2 



Switching between (51) and (49) occurs based on the analysis of exceeding the bound- 
aries Umax of control actions according to the algorithm described in the end of section 
4. 

Informally, the essence of the proposed control is reduced to that the control object 



(41) in area close to 



is parametrically replaced to (47). For this parametrized plant 



conventional feedback linearization is applied. 



5.3 Three phase induction motor control: simulation 

This numerical experiment is conducted to simulate the start of induction motor and 
stabilization of speed and flux. In area 0^ ~ induction machine cannot be linearized by 
feedback and open loop controller usually used for start [18j. With our parametrization 
its possibe to perform feedback lineraization of motor on start as well. 

The model in Simulink to control the system shown in figures 5-6. 
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Figure 5. Siniulink model of overall system. 
The model consists of the following hierarchical blocks: 



AC motor model - model of induction motor in the form (5.3), 



• Power source - inverse Park transformation and the SVPWM signal generator, 

• Controller - subsystem with the PI controllers and the linearizing transformations 
for speed control u flux 0^, 

• Observer - subsystem block implements the calculation of the slip speed UgHp and 
flux (pj. from measured currents z^d, isq 

The controller subsystem is shown in Figure 5. Model of the controller consists of the 
following blocks 



Plant [w, phi] - nonlinear coordinate transformation in form of (49) to control the 
speed uj and 0^, 



Plant [i_s_dq] - nonlinear coordinate transformation in form of (44) to control the 
stator currents isd and isq-, 

PID Controller [i_sd] - Pl-controller for stabilization of isd-, 

PID Controller [i_sq] - Pl-controUer for stabilization of ig^, 

PID Controller [phi_r] - Pl-controUer for stabilization of <^^, 

PID Controller [omega] - Pl-controUer for stabilization of w. 




Figure 6. Controller subsystem model. 
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To increase the adequacy of the modehng to the signals of the stator current isd and 
isq added additive Gaussian noise with variance 0.005 obtained from a source of pseudo- 
random numbers, which corresponds to the 15 mA random error of current measurement. 

Parameters of the model motor are presented in Table 1. 



Symbol 



P 

Rs 
Rj- 

J 



P 

T 

J. rr 



Value 



4 
0.175 

1.2 
0.873 
0.195 
0.195 
0.013 

2 

2 



Description 



rated power, kW 
mutual inductance, H 
stator resistance, ohm 
rotor resistance, ohm 
stator inductance, H 
rotor inductance, H 



m 



combined inertia, kg 

pole pairs 

load torque, N m 

Modeling results are shown on figures 7-9. From the data obtained it is clear that with 
the presence of noise in the current feedback control algorithm provides acceptable perfor- 
mance. In particular, the accuracy of speed control is 0,2 % and accuracy of maintaining 
the magnetic flux is 2 %. 
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Figure 7. Time plot of the motor speed uj{t). The objective of control is smooth 
acceleration to a speed 100 rad/sec within 2 seconds. 
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Figure 8. Time plot of the rotor flux (prit)- The objective of control is maintain the 
magnetic flux of the rotor on the value of 0.31 Wb. 
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Figure 9. Control inputs: the stator voltages Ugd and Usq- 
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Figure 10. Time plots of the stator currents isd{t) and isq{t)- 



6 Three phase induction motor control: experimental 
implementation 

In this section we describe the results of experimental studies to verify and test the 
proposed approach drive control. General view of setup is shown in Figure 11. Three- 
phase asynchronous motor is mechanically connected to the controlled synchronous motor, 
which is used as a torque source. 

As a platform for implementing control algorithms used by the controller dSPACE 
DS5202, which is a hardware target for automatic code generation from MATLAB Simulink 
model. With technology of automatic code generation for hardware targe all implemented 
in Simulink algorithms were tested on a real induction motor without modifications. For 
the motor power supply power-stage based on frequency converter SEW MoviTrac is used. 
PWM control signals for transistors generated in dSPACE controller. 
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Figure 11. Overall view of the experimental setup. 

The experiment was replicated conditions similar to those used in the simulation, 
namely, the braking torque T^ = 4.28 N m. Tracking of speed setpoint ujref{t) was tested 
during experiments. The experimental results of speed control are shown in Figures 8-10. 
Torque was applied at moment t = 32 sec. 
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Figure 12. Time plots of actual speed uj{t) and reference speed Uref{t)- 
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Figure 13. Control inputs: the stator voltages Usd and Ugq- 
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Figure 14. Control outputs after the linearization v^^ and v, 
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Figure 14. Time plots of the stator currents isd{t) and isq{t). 

From these results it is evident that the speed of rotation uj{t) is strictly corresponds 
to the set point Uref{t), the deviation does not exceed 0.8 %. 

7 Conclusion 

In this paper we propose a new method for affine control systems, which combines the 
conceptual simplicity of feedback linearization methods and at the same time expands the 
scope of their applicability to irregular system with poorly expressed relative degree. 

The method tested on an abstract system MIMO sistem with singularities in state 
space. Application of proposed method demonstrated on the speed and rotor magnetic 
flux control in the three-phase asynchronous motor. It has been modeled taking into ac- 
count the effect of measurement errors of the stator currents. Also, the proposed approach 
is implemented and tested on an experimental setup. 

Future work will focus on the investigation of uncertainties influence in an explicit 
form, the generalization of the approach using methods of differential geometry. 
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